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The spinless Salpeter equation supplemented by an instanton induced force is used to describe the 
spectra of light mesons, including the pseudoscalar ones. The coupling constants of the instanton 
induced potential, as well as the quark constituent masses are not treated as simple free parameters 
but are calculated from the underlying instanton theory. Quite good results are obtained provided 
the quark are considered as effective degrees of freedom with a finite size. A further test of the 
model is performed by calculating the electromagnetic mass differences between S-wave mesons. 
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T^j- ■ I. INTRODUCTION 

o . 

Numerous papers has been devoted to the study of meson spectra in the framework of the potential model. In 
most of these works, it is assumed that the quark interaction is dominated by a linear confinement potential, and that 
O . a supplementary short-range potential stems from the one-gluon exchange mechanism (see for instance Refs. 0,|2|)- 
The results obtained with these models are generally in good agreement with the experiments but the mesons r\ and 
rj cannot be described without adding an appropriate flavor mixing procedure with supplementary parameters. 

On the other hand, Blask et al. [3| have developed a non-relativistic quark model which describes quite well all 
mesons (including 77 and rj 1 ) and baryons composed of u, d or s quarks. The long-range part of their interaction is 
the usual linear confinement potential, but their short-range part is a pairing force stemming from instanton effects. 
This force presents the peculiarities to act only on quark- antiquark states with zero spin and zero angular momentum, 
and to generate constituent masses for the light quarks. The main problem of this model, and more generally of all 
non-relativistic models, is that the velocity of a light quark inside a meson is not small compared with the speed of 
light. This makes the interpretation of the parameters of such models questionable 2, p. 164]. 

At present, several works have been devoted to the study of mesons with the instanton induced forces in the 
framework of relativistic or semi-relativistic models but in all the cases constituent masses and coupling 

constants of the instanton induced forces have been considered has free parameters fitted to reproduce at best meson 
spectra. Actually, these quantities can be calculated from instanton theory. In this work, our purpose is to develop 
■ a semi-relativistic model for meson spectra including the instanton induced forces, but with parameters calculated, 
1 as far as possible, with the underlying theory. We will show that such a procedure is possible and gives very good 
JL 1 results provided the quarks are considered as effective degrees of freedom with a finite size. 

Our paper is organized as follows. In Sec. |HJ the model is constructed and the various parameters are presented, 
including the quark size parameters. The numerical techniques and the fitting procedure for the parameters are 
described in Sec. II I II where the results are analyzed and a further test of the model is achieved with the calculation 
of the electromagnetic mass differences. Concluding remarks are given in Sec. II VI 
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II. MODEL 



A. Spinless Salpeter equation 



Our models rely on the spinless Salpeter equation. This equation is not a covariant one, but it takes into account 
the relativistic kinematics. It can be deduced from the Bethe-Salpeter equation by neglecting the retardation effects, 
the mixing with negative energy states and the spinor structure of the eigenstates. The spinless Salpeter equation 
has been often used to describe meson spectra (see for instance Refs. [HljQI)- This equation has the following form 

H = ytp 2 + mf + \jp 2 + mj + V(f), (1) 

where V is the potential between the particles and where p is their relative momentum. The vector p is the conjugate 
variable of the inter-distance f. 

As usual, we assume that the isospin symmetry is not broken, that is to say that the u and d quarks have the same 
mass. In the following, these two quarks will be named by the symbol n (for normal or non-strange quark). 
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B. Funnel potential 

It is now well accepted that the long range part of the interquark interaction is dominated by the confinement. 
The best way to simulate this phenomenon in a semi-relativistic equation is to use a potential increasing linearly 
with the distance. With such an interaction, the Regge trajectories of light mesons are well reproduced 0, P- 137]. 
Moreover, lattice calculations also find that the confinement is roughly proportional to r = \r |. The short range part 
is very often a Coulomb-like interaction stemming from the one-gluon exchange process. The idea is that, once the 
confinement is taken into account, other contributions to the potential energy can be treated as residual interactions. 
It is worth noting that the contribution of a constant potential is always necessary to obtain good spectra. Finally, 
the central potential considered in our models is the so-called funnel potential 

V(r) = -- + ar + C. (2) 
r 

Despite its simple form, this potential was the first one which is able to reproduce the charmonium spectrum quite 
well 0- Moreover, it gives very good results even in the light meson sector (see for instance Ref. 0). 

Experiment shows that vibrational and orbital excitations give much more contributions to the meson masses than 
variations of spin 5 or total angular momentum J in a multiplet. In this case, a spinless Salpeter equation, only 
supplemented by the funnel potential, can yield very satisfactory results. Nevertheless, the situation is completely 
different for the L = mesons for which the mass differences between 5 = and 5 = 1 states is extremely large. For 
these mesons, another interaction must be taken into account in the model. 



C. Instanton interaction 



The instanton induced interaction provides a suitable formalism to reproduce well the pseudoscalar spectrum. 
Indeed it is possible to explain the masses of the pion and the kaon and to describe states with flavor mixing as r\- 
and ?/ -mesons. The form of the interaction depends on the quantum numbers of the state jjjj. 



For L ^ or 5 ^ 0: 



• For L = S = and 1 = 1: 



For L = 5 = and I = 1/2: 



V lnst = 0; (3) 
V lnst = -8g6(f); (4) 
Vi nst = -8g'5(f); (5) 



• For L = S = and I = 0: 

V2g' 

in the flavor space (l/-\/2(|tm) + \dd)),\ss)). The parameters g and g' are two dimensioned coupling constants 
defined as 

3 

.9 = gffcff(s), (7) 
3 

g = g9e&(n), (8) 
ffcffW = Q^ 2 ) dpd (p) P 2 (m^p 2 c l ), (9) 

where ra- is the current mass of the flavor i and a = (2/3)7r 2 (giCft), (qiqi) being the quark condensate for this flavor. 
The function do(p) is the instanton density as a function of the instanton size p. For three colors and three flavors 
this quantity is given by |4|, 

Mp> = 3MW ~i,m)' exp (rm) (10) 
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where 



8tt 2 \ , / 1 

7(p)J =91n U 



9 V Up 



(ii) 



within two-loop accuracy Q . The quantity A is the QCD scale parameter and p c is the maximum size of the instanton. 
This is a cutoff value for which the lnln-term in Eq. is still reasonably small compared with the ln-term. 
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contributions to the constituent masses. The expression of these contributions are given by [j 



An interesting property of the instanton induced interaction is the renormalization of quark masses, as it gives 

T 

4 

Am,, = — 
3 

and 



dp d (p) (m° n - p 2 Cn) (m° - p 2 c s ), 



(12) 



4 [ Pc 

-it 2 J dp d {p) (m„ - p 2 c n f 



(13) 



The instanton interaction is not necessarily the only source for the constituent masses 3]. Actually we introduce two 
supplementary terms S n and S s which can be added to the running masses. These terms are free parameters and 
are not dependent on the instanton parameters. Results with vanishing and non- vanishing S n and/or S s are given in 
Table ITT1 Finally the constituent masses in our models are given by 



m„ 

m. 



(14) 
(15) 



We can rewrite expressions |JS}, (|12ll and (|13|l in a more interesting form for numerical calculations by setting a 
dimensionless instanton size 



and defining another dimensionless quantity 



where x c = Ap c . So we obtain 

9 = 
9' = 
Am„ = 



dx 



Sir 2 1 
2 A 3 
Sir 2 1 



9 In [ - 
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x = Ap, 



In [ - 

x 



-32/9 



2 A 3 

.1 
•a 



m° s au(x c ) - a 13 (x c ) 
m° an(xc) - ai 3 (x c ) 



m„m" a 9 (x c ) 



(c n m° s + c s m°) 
A 2 



an(x c ) + 



c n c s 



ai3{x c ) 
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^ °) 2 a 9 (x c )-2 



( \ 2 

n an(x c ) + ai 3 (x c ) 



A 2 



(16) 
(17) 

(18) 
(19) 
(20) 

(21) 



with S = 3.63 10~ 3 x 47r 2 /3. Except the quantity x c , all parameters involved in Eqs. I|18 |l -I|21 |l have expected values 
from theoretical and/or experimental considerations. The integration in Eq. 1)17(1 must be carried out until the ratio 
of the ln-term on the lnln-term in Eq. (|llfl stays small . This ratio, called R here, increases with x from zero at 
x = x\ = 1/e to very large values (see Fig. 01 . At x = X2 ~ 0.683105, the value of this ration is 1. This last value 
corresponds to the minimum of the instanton density (see Fig. 0) . Thus we define the parameter e by 



x c = X\ + e(x2 — xi) with e 6 [0, 1]. 



(22) 



In this work e is a pure phenomenological parameter whose value must be contained between and 1. To save 
calculation times, we have calculated some values of a n (x c ) functions for some given values of e, and we use a spline 
algorithm to find others values. A good accuracy can be obtained since these functions and their derivatives are 
known. We have verified that the spline procedure allows an accuracy better than 10 -3 on values of quark masses 
and instanton coupling constants. So, no numerical integration is necessary. The functions ag(x), an(x) and 0:13(2;) 
between xi and X2 are given in Fig. [21 
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D. Effective quarks 



The quark masses used in our model are the constituent masses and not the current ones. It is then natural to 
suppose that a quark is not a pure point-like particle, but an effective degree of freedom which is dressed by the 
gluon and quark-antiquark pair cloud. As a correct description of this effect is far from being obvious, we use a 
phenomenological Ansatze, as it is the case in many other works (see for instance Refs. It seems natural to 

consider that the probability density of a quark in the configuration space is a peaked function around its average 
position. The form that we retain is a Gaussian function 

MO = ^3/2 eM-rVji). (23) 

It is generally assumed that the quark size ji depends on the flavor. So, we consider two size parameters 7„ and 7 S 
for quarks n and s respectively. Any operator which depends on the quark positions f and fj must be replaced by an 
effective one which is obtained by a double convolution of the original bare operator with the density functions pi and 
Pj. As a double convolution is a heavy procedure which generates very complicated form for the convoluted potentials, 
we assume that the dressed expression Oijif ) of a bare operator Oij(f), which depends only on the relative distance 
f — fi — fj between the quarks qi and qj , is given by 



On (r ) = J df {j (f') Pij (f -f), (24) 
where pij is also a Gaussian function of type l|23[l with the size parameter 7^ given by 



7*i = V^+%' ( 25 ) 

This formula is chosen because the convolution of two Gaussian functions, with size parameters 7^ and jj respectively, 
is also a Gaussian function with a size parameter given by Eq. (JS5J. 

After convolution with the quark density, the funnel dressed potential has the following form 



2 



lij eM~r 2 /%) ( | 7 tJ 
JUr \ 2r 2 



erf(r/7y) 



C, (26) 



while the Dirac-distribution is transformed into a Gaussian function 



6(f) = exp(-r 2 /7y-)- (27) 

Despite this convolution, we will consider, for simplicity, that the instanton induced forces act only on L = states. 

The modification of the confinement potential seems very important but, actually, only its short range part is 
modified. It has little effect for a potential whose essential role is to govern the long range dynamics. It has been 
verified @ that the convolution of the linear potential could be neglected without changing sensibly the results, but 
we have nevertheless used the form <|26|) to be consistent. 

The Coulomb part of the interaction is transformed into a potential with a finite value at origin. This can be seen 
as a mean to simulate asymptotic freedom. It has also been noticed that the use of an error function remove the 
singularity of the spin-spin interaction, when this correction is taken into account 0, p. 162]. 

The introduction of a quark size is only necessary, from a mathematical point of view, to avoid collapse of the 
eigenvalues due to the presence of a Dirac-distribution in the instanton induced interaction. The definition of the 
size parameter 7^ is not obvious in the case of the non-diagonal term of the instanton induced interaction ©. This 
matrix element mixes \nn) and \ss) flavor states. So, what definition choose for 7^ , that we will note 7 m f in this case? 
A first possibility is to take 



7mf = VlZ+ll (28) 
as in the case of \sn) mesons. But we can also choose 

7mf = V^lnls, (29) 



for instance. Other choices are possible, but we only take into account these two definitions in the following. 
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III. NUMERICAL RESULTS 
A. Numerical techniques 

The eigenvalues and the eigenvectors of our Hamiltonian are obtained in expanding trial states in a harmonic 
oscillator basis \nlm). In such a basis, matrix elements of the potential are given by 

Z+n+n 

(n'lm\V(r)\nlm) = B(n' ,l,n,l,p)I p , (30) 

P =i 

with 

2 f°° 

Ip = T(p + 3/2) J dxx2p+2e M-x 2 )V(bx). (31) 

The quantities I p are the Talmi integrals, while the coefficients B(nf ,l,n,l,p) are geometric factors 8] which can be 
calculated once for all. The parameter b is the oscillator length which fixes the scale of the basis states. With this 
parameter a dimensionless length x — b/r can be defined. In our model all the Talmi integrals for the potential part of 
the Hamiltonian are given by analytical expressions |(| . The matrix elements of the kinetics part can be calculated with 
numerical integrations, but we prefer to use another technique, much accurate and less time consuming. This method 
is described in Ref. 0. We just give here the main hits of the procedure. We calculate the matrix elements of the 
operator p 2 + m 2 in the oscillator basis (this is an analytical calculation) . Let us assume that P is the corresponding 
matrix, D the diagonal matrix of the eigenvalues of P, and U the transformation matrix (D = U^ 1 P 77). Then 
the matrix elements of the operator \/p 2 + m 2 , are contained in the matrix U D 1 / 2 U~ l . This procedure is about 
six times faster than a numerical integration and gives a better accuracy. The utilisation of this procedure and the 
analiticity of the Talmi integrals lead to very short time of calculation. A minimization with thirteen parameters lasts 
about 10-20 minutes on a Pentium 200 workstation. 

At last, it is worth noting that all the results obtained with the technique described above have been verified with 
the three-dimensional Fourier grid Hamiltonian method 0. 

B. Fitting procedure 

The purpose of this work is to try to reproduce the spectrum of light mesons with a quite simple model. Indeed 
we use only a central potential supplemented by an instanton induced interaction to describe the pseudoscalar sector. 
We need thirteen parameters to obtain a satisfactory theoretical spectrum. The instanton interaction is defined by six 
parameters: The current quark masses, the quark condensates for the flavors n and s, the QCD scale parameter A, 
and the maximum size p c of the instanton. Three parameters are used for the spin-independent part of the potential: 
The slope of the confinement a, the strength k for the Coulomb-like part, and the constant C which renormalize the 
energy. Four supplementary parameters are introduced: The effective size of the quarks n and s, and two terms 5 n 
and S s , which contribute to the constituent quarks masses. 

In our model, the quantum numbers L, 5 and I are good quantum numbers. For L = states, the instanton induced 
interaction raises the degeneracy between 5 = (pseudoscalar) and 5 = 1 (vector) mesons, but 5 = 1 nn-mesons, 
which differ only by isospin, have the same mass. This is in quite good agreement with the data. In the sector L ^ 0, 
instanton induced interaction vanishes and the potential is spin independent. Consequently, states which differ only 
by 5, J and I are degenerate. In first approximation, this corresponds also quite well to the experimental situation. 

To find the value of the parameters, we minimize a x 2 function based on the masses of 18 centers of gravity (c.o.g.) 
of multiplets containing well-known mesons (see Table [Q: 

i 

where the quantity AM™ P is the error on the experimental masses (it is fixed at the minimum value of 10 MeV). 
The quantum number v (vibrational or radial quantum number), L, and 5 of these mesons are determined on an 
assignment made in the Particle Data Group tables [l(J]. It is worth noting that all members of several multiplets 
considered here are not known. So we calculate the center of gravity with only the known mesons; no attempt is made 



M 



M oxp 



AAl cxp 



(32) 
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to estimate the masses of the missing states. The usual procedure to define a center of gravity is 

Mc -°^- E^(2/+1)(2J+1) ' 

In the following, a center of gravity will be indicated by the name of the state of the multiplet with the higher quantum 
numbers J and /. At last, note that we do not include the mesons ao(980) and /o(980) in the L = 1 multiplet since 
experimental considerations and some theoretical works [iH p. 99, 557] suggest that they are not gg-mesons. 
To perform the minimization, we use the most recent version of the MINUIT code from the CERN library ^ 



C. Meson spectra 



A great number of parameter sets have been found for our model Hamiltonian. We present here only five among 
the most interesting ones, denoted from I to V. Parameters for these models are presented in Table |n] All these 
models have common features. For instance, it is always possible to obtain a quite good fit with a low current quark 
mass m°, in agreement with the bounds expected 0. The strange quark mass m° s varies more significantly following 
the model, but its value is always reasonable ^(j- The QCD scale parameter A is a very sensible parameter of our 
model since a small change can largely modify the quantities deduced from the instanton theory (see Eq. (|18|l to J2J). 
Nevertheless, the values found are always in agreement with usual estimations [Tof. Values for the quark condensates 
are also reasonable |12|. There is no direct measurement of the string tension parameter a, but lattice calculations 
favor the value 0.20 GeV 2 with about 30 % error [13j. Our values are in good agreement with this prediction. 

The strength k of the Coulomb-like potential can strongly vary from one model to another (see for instance Ref. 0). 
The values found in this paper are always less than 0.5, which can be considered as an upper limit for relevant values. 
In model V, we have fixed k = as in the model of Ref. But the \ 2 obtained is not good; this is essentially due 
to radial excitations: p(1450), 7r(1300) and 0(1680). We conclude that the Coulomb-like potential must be taken into 
account and that instanton induced forces cannot explain alone the short range part of the potential. Actually, this 
remark is already mentioned in the Ref. Q. As we can see from Tabled the constant potential is always necessary to 
obtain good spectra. Its origin is not simple to explain. It can be considered as a mechanism linked to the flux-tube 
model of mesons [f|. It has also been suggested that C « 2, p. 190]. Our models are not in very good 

agreement with this last prediction, as we can note deviation as large as 40 %. But we do not consider that it is as 
an important drawback of the models. 

There is no real estimation of the quark sizes. They are purely phenomenological parameters whose role is to take 
into account relativistic effects 0] as well as sea-quark contributions. In all our models the size of the n-quarks 7„ 
is nearly a constant, around 0.7-0.8 GeV -1 . The n-quark size drops down to 0.6 GeV -1 in model V without the 
Coulomb-like potential. This shows the strong influence of this interaction in determining 7„. It is worth noting 
that the introduction of quark sizes modifies deeply the structure of the potential. This is illustrated on the Fig. 
where the dressed and non-dressed potential for the p-meson are presented. The s-quark size depends on the Ansatze 
chosen to calculate 7 m f. When 7 m f = a/tJ+tJ, 7s can take vanishing values without generating collapse of the 
eigenenergies. This is no longer true if 7 m f = ^/2j n -f s . The first definition of 7 m f is chosen in model II and we can 
remark that in this case the value of 7 S is significantly smaller than in other models, where the second definition is 
used. 

Table [2] shows that the contributions of parameters 5 n and 8 S can be quite large when they are not fixed in the 
minimization. The origin of these terms is not clear but their values indicate that the instanton effects cannot generate 
solely the constituent masses in our model Hamiltonian. We have found sets of parameters which give reasonable 
X 2 values but, in this case, the constituent quark masses can be very small. Model III has actually the lowest \ 2 
value that we have found but the constituent quark masses are so small that a generalization of this model to baryon 
spectra appears very problematic |l4j . 

In all sets of parameters that we have determined, the value of e is always close to zero. This is consistent with the 
fact that the cutoff radius for the integration over instanton density must be small enough in order that the lnln-term 
in Eq. (|llfl must be reasonably small compared with the ln-term. Actually, we can fix arbitrarily e = without 
spoiling the results, while values close to unity yield bad results. 

In Figs. 0] to El we present the meson spectra of model I. We consider this model as the best we have obtained. 
Model III is characterized by a lowest % 2 value but, as mentioned above, the corresponding constituent quark masses 
are too small to hope to obtain good baryon spectra. The \ 2 value of model II is also good but the strange quark size 
seems very small with respect to usual estimations found in the literature. In order to decide between models I and 
II, we have tested these two potentials in the heavy meson sector. By fitting the masses and the sizes of the quarks c 
and b on 10 c.o.g. of cc, bb, sn, cs and bn mesons, we have recalculated a \ 2 value for the all 28 c.o.g. considered in 
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this paper (18 for light mesons and 10 for heavy mesons). The results are indicated in Table ITU Clearly, model I is 
preferable. The bad value obtained for model II is due to an interplay between the small strange quark size and the 
low Coulomb strength. 

The partial \ 2 f° r each meson (actually, each c.o.g. of a meson multiplet) of model I is below 1.6, except for the 
7r(1300)-, p(1450)- and a4(2040)-mesons. In Fig. we can see that the nn mesons are quite well reproduced. The 
larger error is obtained for the 7r(1300), but the uncertainty on this meson is very large. We can remark slight 
deviation from the linear Regge behavior, but it can be very well reproduced within semi-relativistic kinematics when 
only a linear potential is used. 

In the sn sector, vibrational excitations of the L = mesons are not satisfactorily reproduced (these states are not 
taken into account in the minimization procedure). From the experimental point of view, the situation is not clear 
[To). For instance, internal quantum numbers of if* (1410) and of if* (1680) are not well defined yet. This problem 
was also revealed in some previous works [H [f| . But all the ground states are quite well reproduced. 

The ss sector is poor in experimental data but all the states obtained in our model are in quite good agreement 
with these data. In the mixed flavor mesons, we reproduce the rj and rf states. The two vibrational excitations of 
these mesons can be identified quite well with the r/(1295) and the 7/(1760). Note that a calculated state lies between 
the 77(1295) and 7/(1440), which is a non-qq meson candidate [TJ}. At last, we find a supplementary state between the 
77(1760) and the 7/(2225), the last one being a not-well established state. 



D. Electromagnetic corrections 



The electromagnetic mass differences between mesons are usually ignored in studies of meson spectra. The reason 
is that these mass differences are very small, of the order of some MeV, compared with the orbital and vibrational 
excitations due to the strong interaction, which can amount about one GeV. Nevertheless, it is interesting to calculate 
the electromagnetic splittings as they provide further tests on the model. In particular, the mass differences are very 
sensitive to the short range part of the wave functions. 

The electromagnetic mass differences between mesons are due to two distinct effects. A first contribution is provided 
by the mass difference between the u and d quarks. In our models we have assumed that these two quarks have the 
same mass, but it is no longer relevant to calculate electromagnetic phenomena. In the following, we will assume that 
m n = (fhd + fh u )/2, where fhi is the real constituent quark mass. We define ed — fhd — m n and e u = fh u — m n . For 
the strange quark, we have fh s = m s and e s = 0. The second contribution is due to the electromagnetic potential 
existing between quarks. In first order approximation, this interaction has the following form 



V cm (r) 



QiQj 



(34) 



where a is the electromagnetic fine-structure constant, and Qi is the quark charge for the flavor i in unit of e. This 
approximation is too crude since the first relativistic corrections of V em and V em have similar contributions to the 
electromagnetic mass differences 0. We will consider only S-wave mesons, that is to say that the spin-spin interaction 
is the only non- vanishing correction. As we work in the framework of semi-relativistic models, we will use a relativized 
version of the spin-spin potential. A method to obtain such a potential is suggested in Ref. 0, P- 198] and has been 
used in Ref. [l] . The idea is to replace the factors 1 /in; by the operators 1 / \f p 2 + mf in the interaction expression. 

The contribution of the electromagnetic Hamiltonian H em can be calculated within the first order in perturbation 
theory • Assuming that this Hamiltonian is the difference between the total Hamiltonian (strong plus electromag- 
netic) and the Hamiltonian used in our models, the electromagnetic contribution is in first order approximation 



(H r . 



1 



+ m,- e 



3 *-o 



1 



a— ( QiQj Si ■ S 3 



-Mr 



(35) 



where the symbol [O] is used to indicate the symmetrization of the operator O. The expression (|35|) reduces to the 
usual non-relativistic form when quark masses tend toward infinity 0, P- 169]. Mean values of operators can be 
evaluated as well by using the harmonic oscillator development of the wave functions as well directly in configuration 
space. The isospin breaking between the u and d quarks is a free parameter in our models, so this quantity has been 
fitted to reproduce the K + -K° mass difference. The results for model I are shown in table ITTTI Other models give 
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similar results since in all cases the pseudoscalar states are well described. In magnitude, all results are around two 
times the experimental data, but the hierarchy of mass differences is preserved. In particular the quantity m c i — m u 
has the good sign. This is not always obvious to obtain in potential models [lfij |. It is worth noting that our models 
are optimized to reproduce meson mass spectra with an accuracy of around 10 MeV, and not to obtain the best 
possible wave function. So, we consider that the results found are quite well satisfactory. 

We also tried to calculate the electromagnetic mass differences with the non-relativistic equivalent of Eq. I[ri5[) . but 
the results obtained are very bad. For instance, isospin breaking between d and u quark can be found as large as 
50 MeV. We have found that the non-relativistic form of the electromagnetic spin-spin interaction is responsible of 
so poor results. This indicates that relativized forms of operators are important to obtain coherent results in semi- 
relativistic models. Since the Coulomb-like strong interaction is of pure vector type, a strong spin-spin interaction 
must be introduced if one wants to take into account the relativistic corrections of the strong potential. We can guess 
that the use of a relativized form will be highly preferable than the usual non-relativistic form. 

IV. CONCLUDING REMARKS 

It is well known that the spinless Salpeter equation supplemented by the so-called funnel potential can describe 
quite well the main features of the meson spectra Q- The mass differences between members of a L-multiplet are 
generally small with respect to vibrational and radial excitations. Except in the pseudoscalar meson sector, the 
spinless approximation is then generally good. The instanton induced forces provide a satisfactory way to describe 
the structure of light S-wave mesons, including the annihilation phenomenon. 

Several works are devoted to the study of meson spectra within relativistic or semi-relativistic models using instanton 
effects In these works, the instanton parameters are simply considered as free parameters to be fitted on 

data. In our paper, we try to calculate these parameters from the underlying theory. Our model relies on the 
spinless Salpeter equation supplemented by the usual funnel potential and an instanton induced force. Thirteen 
parameters are necessary to completely fix the Hamiltonian. This number could appear large, but six parameters 
are strongly constrained by experimental or theoretical considerations, namely, the current quark masses, the QCD 
scale parameter, the quark condensates and the string tension. The Coulomb strength and the constant potential are 
unavoidable parameters of potential models. Actually, the Coulomb-like potential could be replaced by an interaction 
taking into account the asymptotic freedom, but it is a complication especially necessary when describing heavy meson 
spectra. 

The relevance of the remaining five parameters is most questionable. If quarks are considered as effective degrees 
of freedom, then a size as well as a constituent mass can be associated with a quark. Generally, the quark size is 
described with a two parameter function of the constituent mass. As wc consider only two different flavors, it is useless 
to introduce such a quark size parameterization. Up to our knowledge, there is no reliable estimation of the quark 
size, but our values are in good agreement with values found in other works 0, ||| . The radial part of the instanton 
induced potential is a Dirac-distribution, so it is necessary to replace this form by a short range potential peaked at 
origin. The introduction of quark sizes offers a natural way to realize that, without introducing a new free parameter 
which is the instanton range. 

The two dimensioned coupling constants of the instanton induced force are calculated from the instanton theory. 
This interaction generates also constituent masses for the light quarks. Our work shows that other contributions must 
be added to explain the large constituent masses of n and s quarks. These contributions which lay between 100 and 
200 MeV are quite large. It is possible to cancel them, but the price to pay is then to obtain unrealistic very small 
constituent masses. We consider that it is better to keep the large values of these contributions. Their origin is not 
clear but sea-quark or relativistic effects could explain their presence. 

The parameter e reflects the uncertainties about the form of the instanton density, which is known within the 
two-loop approximation. The values found in all models are consistent with the fact that the cutoff radius of the 
instanton size must be small. 

We have shown that the light meson spectra can be described by a semi-relativistic model including an instanton 
induced force whose characteristics are calculated from the underlying theory. Moreover, calculation of electromagnetic 
mass differences between mesons indicate that the wave function obtained are quite satisfactory. The next steps are 
to calculate meson widths and to generalize the model to baryons. Such a work is in progress. 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

x = Ap 

FIG. 1: One- loop and two- loop approximations of the instanton density (see Eq. (IIUI '1 as a function of the dimensionless 
instanton size x — Ap. The absolute value of the ratio R (see Sec. Ill Cl . and the values of x corresponding to e = and e = 1 
are indicated. 
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TABLE I: Centers of gravity (c.o.g.) of L- and I-multiplets in GeV for mesons chosen to fix the parameters of the models. 
The values of the c.o.g. and their corresponding errors are given by formula J33|> . The symbol "mf " means "mixed flavor". A 
meson name used to represent a multiplet in Figs. [I] and [(J is underlined. 

State Flavor I J p(c) 7V 2S+1 Lj c.o.g. 

7T nn 1 (T+ 1 0.138±0.003 



UJ 


tin 





1~ 


1 


3 S 1 


0.772±0.001 


p 


nn 


1 


1~ 


1 


3 Si 




fti(1170) 


nn 





1+- 


1 


'Pi 


1.265±0.013 


6i(1235) 


nn 


1 


1+- 


1 


l Px 




/i(1285) 


nn 





1 ++ 


1 


Pi 




ai(1260) 


nn 


1 


1++ 


1 


3 Pl 




/a (1270) 


nn 





2++ 


1 


3 P2 




a2(1320) 


nn 


1 


2 


1 


3 D 
P2 




7^(1670) 


nn 


1 


2 


1 


1 n 
L>2 


l.o81±U.012 


w(1600) 


Tin 





1~ 


1 


3 D, 




p(1700) 


nn 


1 


1" 


1 


3 D 1 






1 VI l 


n 
u 


q 


i 
i 


^3 




p 3 (1690) 


nn 


1 


3~ 


1 


3 D 3 




/ 4 (2050) 


nn 





4++ 


1 


3 F, 


2.039±0.022 




mi 


1 
1 


A++ 


1 

1 


3 F< 

FA 




tt(1300) 


nn 


1 


0"+ 


2 


'S 


1.300±0.100 


cj(1420) 


nn 





1~ 


2 


3 S 1 


1.454±0.026 


p(1450) 


nn 


1 


1~~ 


2 


3 




K 


sn 


1/2 


0" 


1 




0.496±0.002 


K* (892) 


sn 


1/2 


1" 


1 


3 S 1 


0.892±0.001 


Ai(1270) 


sn 


1/2 


1+ 


1 


'Pi 


1.382±0.005 


A ( *(1430) 


sn 


1/2 


h 


1 


3 r> 

Po 




^(1400) 


sn 


1/2 


1+ 


1 


3 Pl 




A 2 *(1430) 


sn 


1/2 


2+ 


1 


3 P2 




A 2 (1770) 


sn 


1/2 


2~ 


1 


1 n 
P>2 


1.774±0.012 


A'* (1680) 


sn 


1/2 


1" 


1 


3 D 1 






sn 


1/2 


Z 


1 


3 n 

L>2 




A 3 *(1780) 


sn 


1/2 


3" 


1 


3 D 3 




t 


SS 





1"- 


1 


3 S 1 


1.019±0.001 


hi (1380) 


ss 





1+- 


1 


'Pi 


1.482±0.009 


/i(1510) 


ss 





1++ 


1 


3 Pi 




/a (1525) 


ss 





2++ 


1 


3 P2 




03(1850) 


ss 





3"- 


1 


3 D 3 


1.854±0.007 


0(1680) 


ss 





1~~ 


2 




1.680±0.020 


/a (2010) 


ss 





2++ 


2 


3 P2 


2.011±0.080 



-f n n — h ilo n 1 n nm 



11 



TABLE II: Optimal parameters found for various models (see Sec. 1111 CI . Values between square brackets are fixed quantities 
in the model. When available, the expected value of a parameter is also given in the column "Exp.". The Ansatze to calculate 
7 m f (see Sec. Ill Dl is indicated and the parameters m n , rn Sl g, g and p c are calculated. The quantity x 2 (light) is the value of 
the x 2 given by relation (13211 for the set of mesons from Table U] The corresponding quantity for an extended set of mesons 
(see Sec. 1111 d is indicated by x 2 (ah)- 



Parameter 


Model I 


Model II 


Model III 


Model IV 


Model V 


Exp. 


m° (GeV) 


0.015 


0.013 


0.004 


0.015 


0.015 


0.002-0.015[10] 


m° (GeV) 


0.215 


0.196 


0.199 


0.271 


0.166 


0.100-0.300[10] 


A (GeV) 


0.245 


0.299 


0.223 


0.231 


0.262 


0.209in 3 33[10j 


(nn) (GeV 3 ) 


(-0.243) 3 


(-0.223) 3 


(-0.236) 3 


(-0.250) 3 


(-0.240) 3 < 


(-0.225 ±0.025) 3 [12J 


{ss) /{fin) 


0.706 


0.834 


0.704 


0.700 


0.703 


0.8 =b 0.1[12J 


a (GeV 2 ) 


0.212 


0.225 


0.207 


0.201 


0.235 


0.20 ± 0.03[13J 


K 


0.440 


0.283 


0.366 


0.367 


[0.000] 




G (GeV) 


-0.666 


-0.781 


-0.593 


-0.551 


-0.860 




7n (GeV 1 ) 


0.736 


0.730 


0.764 


0.773 


0.586 




7s (GeV" 1 ) 


0.515 


0.186 


0.552 


0.524 


0.409 




S n (GeV) 


0.120 


0.147 


[0.000] 


[0.000] 


0.120 




Ss (GeV) 


0.173 


0.228 


0.118 


[0.000] 


0.194 




e 


0.031 


0.162 


0.011 


0.046 


0.003 




7mf 


\/27n7 s 


+ 7l 


y/2j„ 7s 




7s 




m„ (GeV) 


0.192 


0.206 


0.061 


0.090 


0.176 




m s (GeV) 


0.420 


0.443 


0.351 


0.312 


0.385 




g (GeV" 2 ) 


2.743 


2.767 


3.117 


3.336 


2.002 




g> (GeV" 2 ) 


1.571 


1.213 


1.846 


1.807 


1.241 




p c (GeV" 1 ) 


1.541 


1.402 


1.664 


1.598 


1.461 




X 2 (light) 


14.9 


16.7 


10.9 


73.1 


43.9 




X 2 (all) 


35 


95.5 











TABLE III: Electromagnetic mass differences for some mesons in MeV. The results of our model I is compared with the data. 



Model I Experiment [10] 


- 


- m^o 


7.184 4.594 ±0.001 


"V - 


- m p0 


1.361 -0.300 ±2.200 


m K + 


— m K o 


fitted -3.995 ±0.034 


m K ,+ 


— m K t,o 


-11.047 -6.700 ±1.200 


m d - 


m u 


26.610 a < 15 



"fitted to reproduce the m K + — m K o difference. 
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FIG. 2: Values of the function a n (x c ) (see Eq. 117H ) as a function of x c and e for the three interesting values of n. The curves 
begin at x c = xi = 1/e (e = 0), and end at x c = X2 (e = 1). 




-5 ^ — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — 

1 2 3 4 5 

r(GeV"') 

FIG. 3: Interquark potential V(r) in GeV as a function of r in GeV -1 for the p- and the 7r-meson. The potential for the 
p-meson without taking into account the effect of the quark sizes (non-dressed) is also presented. The corresponding potential 
for the 7r-meson cannot be shown because the presence of a Dirac-distribution coming from the instanton induced interaction. 
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O Experience 
♦ Modele I 



a 7T 2 (2ioo) 



6 P,(2250) 



»|a 4 (2040) 



9 7i(1800) 



e umo) 



» P,(1690) 



ft I P( 14S0)| 
fj |n(1300)| 9 1^(1320)1 



L=0, S=0 L=0, S=l L=l L=2 L=3 L=4 



FIG. 4: Comparison of experimental (open circles with the corresponding error bars) and calculated (black diamonds) spectra 
of nn-mesons for the model I. Framed names indicate centers of gravity of multiplets used to fix the parameters. 



<D 2 

O 



O Experience 
• Modele I 



® K,(2250) 



Q K 3 (2380) 



Q K(I830) 



Q K, (1980) 



ft K 4 (2045) 



O K (1680) 



S Kg (1780) 



OK(1460) OK (1410) S |K, (1430) 
«|k'(892)| 



L=0, S=0 L=0, S=l L=l 



L=2 



L=3 



FIG. 5: Same as for Fig. 0|but for sn mesons. 
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O Experience 
♦ Modele I 



9 f 4 (2220) 



« r,(2225)? ♦ 

§ |f 2 (2010)| 

m [W 1 850)1 
9 * 176 °> ®g(T68og 

q r|(1440)? • 

|f 2 '(l525)| 

O n(1295) 



L=0, S=0 L=0, S=l L=l L=2 L=3 L=4 



FIG. 6: Same as for Fig.|l]but for ss and mixed-flavor mesons. 



